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Abstract 


We investigate the behaviors of entanglement entropy in the holographical 
insulator/superconductor phase transition. We calculate the holographic en- 
tanglement entropy for two kinds of geometry configurations in a completely 
back-reacted gravitational background describing the insulator/superconductor 
phase transition. The non-monotonic behavior of the entanglement entropy is 
found in this system. In the belt geometry case, there exist four phases char- 
acterized by the chemical potential and belt width. 


1 Introduction 


The AdS/CFT correspondence [1, 2, 3, 4] provides a powerful theoretical method to study 
the strongly coupled systems in various fields of physics. One of the most widely investi- 
gated objects is the holographic superconductor (superfluid) [5, 6, 7, 8, 9]. The physical 
picture is that some gravity background will become unstable as one tunes some parameter, 
such as temperature for black hole and chemical potential for AdS soliton, to develop some 
kind of “hair”. This condensation of the “hair” induces the symmetry breaking, which 
results in the non-vanishing vacuum expectation value of the dual operator in the field 
theory side. 

On the other hand, as a measurement of how a given quantum system is entangled or 
strongly correlated, the entanglement entropy is also considered as a useful tool for keeping 
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track of the degrees of freedom of strongly coupled systems. Dividing a quantum system 
into a subsystem A and its complement, the entanglement entropy of A is known as the 
famous von Neumann entropy 


Sa = —Tra(lĝaln ĝa), (1) 


where py is the reduced density matrix for A by tracing over its complement of the total 
density matrix. However, the calculation of entanglement entropy for a given system is 
found to be very difficult except for the case in 1 + 1 dimensions. 

Inspired by the Bekenstein-Hawking entropy of black hole and motivated by the de- 
velopment of AdS/CFT correspondence, the authors in Refs. [10, 11] have presented a 
proposal to compute the entanglement entropy of conformal field theories (CFTs) from 
the minimal area surface in gravity side. Precisely, in an asymptotically anti de-Sitter 
(AdS) spacetime, consider a slice at constant radial coordinate r = 1/e, which is dual to a 
field theory “living” on this slice with a UV cutoff labeled by the small parameter e. the 
€ — 0 corresponds to the UV limit of the dual field theory, which is known as the UV-IR 
relation [12]. Search for the minimal area surface y4 in the bulk with the same boundary 
OA of a region A on that slice. The entanglement entropy of A with its complement is 
defined as the “area law”: hee) 

rea(ya 
sa = Tw, (2) 
where Gy is the Newton’s constant in the bulk. A proof for this proposal from the basic 
principle of holographic duality is given in Ref. [13] and further properties of the holographic 
entanglement entropy are discussed in Ref. [14]. 

This proposal provides an elegant and executable way to calculate entanglement en- 
tropy of a strongly coupled system which has a gravity dual. Since then a lot of works have 
been carried out [15, 16, 17, 18, 19, 20, 21, 22, 23] for investigating and calculating the 
entanglement entropy in various gravity theories. In particular, Refs. [20, 21, 22] presented 
the calculations of entanglement entropy for AdS soliton geometry and Ref. [23] considered 
the case with higher derivative corrections. Very recently, the authors of Ref. [24] have 
studied the entanglement entropy in the holographic metal/superconductor phase transi- 
tion. It was found that the entanglement entropy in superconducting phase is always less 
than the one in the metal phase. The entanglement entropy is continuous, but does not 
smoothly cross the second order phase transition. 

The holographic insulator /superconductor phase transition was first studied in Ref. [9], 
where the AdS soliton background [25] was used to mimic the insulator /superconductor 
phase transition at zero temperature. Motivated by Ref. [24] we are going to explore the 
behaviors of the entanglement entropy in the insulator/superconductor phase transition. 

The model of the holographic insulator/superconductor phase transition is constructed 
in the Einstein-Maxwell-scalar theory with a negative cosmological constant in five-dimensional 
spacetime. After completely solving the set of coupled equations of motion, we calculate 
the entanglement entropy for two geometry configurations, one is a half space and the 
other is a belt. In both cases, we find a discontinuity in the slope of the entanglement 


entropy as a function of chemical potential across the phase transition. But, unlike the 
case of the metal/superconductor phase transition [24], after condensate, the entropy first 
rises as we increase the chemical potential u, arrives at its maximum at a certain chemical 
potential and then it decreases monotonously. For the belt geometry, there exits a so-called 
“confinement /deconfinement” transition controlled by belt width. Our results confirm the 
prediction in Ref. [21] that this phase transition must appear in any consistent gravity 
dual of a confining theory. So there are totally four “phases” parameterized by chemical 
potential u and belt width /. The complete “phase diagram” is drawn in Figure.(8). Here 
we would like to address that in fact there does not exist the parameter “belt width @” in 
the dual field theory. When one changes £, that a new “phase” occurs in the entanglement 
entropy manifests a length scale of dual field theory under consideration. For example, in 
this paper the length scale is nothing, but the correlation length of the confining theory. 

This paper is organized as follows. In Section (2), we briefly review the holographic insu- 
lor/superconductor phase transition presented in Ref. [26], and give the complete equations 
of motion to be solved. In Section(3), the fully back-reacted system is solved by shooting 
method. In Section(4) we explore the behaviors of the entanglement entropy in insula- 
tor/superconductor phase transition. The conclusions and some discussions are included 
in Section (5). 


2 AdS Soliton Background 


Let us start with the Einstein-Maxwell-complex scalar theory with a negative cosmological 
constant in five-dimensional spacetime 


12 #1 
S= | deya (R+ Ta- F” Ew IDU- msl) (3) 
where L is the radius of AdS spacetime, field strength F = dA and D, = V,,—iqA,. When 
the Maxwell field and scalar field are vanishing, the above action admits an exact solution 
named as AdS soliton 


L?dr? 4 
O +1°(—dt? + dx? + dy? + f(r)dy?), f(r) =1- 2. (4) 


This asymptotically AdS solution approaches to R}? x St near the boundary at r — oo. 
Moreover, the compactification x ~ x + 7L/rp needed to avoid a conical singularity at 
r = ro gives a dual picture of a three-dimensional field theory with a mass gap, which 
resembles an insulator in the condensed matter physics [9]. The geometry in the (r, x) 
sector just looks like a cigar with its tip at r = ro. The temperature associated with the 
soliton background is zero. 

In order to mimic the superconductor phase, we need the soliton solutions with nonva- 
nishing w. Furthermore, the back reaction of the matter fields must be taken into account 
in order to see the difference of the holographic entanglement entropy. So we choose the 


ds? = 


ansatz as follows 


dr? 


2 __ 
oe r?B(r) 


+r’ (-e dt? + dr? + dy? + eA B(r)dx”), Ar=¢(r), W=V(r). (5) 


Here we have set L = 1 without lose of generality. We require that B(r) vanishes at the 
tip of the soliton. And in order to obtain a smooth geometry at the tip ro, x should be 
made with an identification 
Ate Alro)/2 
~x4+T T= —.—_.. 6 
x~x+T, Bro) (6) 


The independent equations of motion under this ansatz are given as 


> 5 A! B! C' e7 242 
w+ (2 +5425 S)v tag 5 (=n) =o, (7) 
3 A! B' C! 2? Q 
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r(6+rC') í 
n la Dg BV. op 2g py? 7 _ 
OP Ct te | Oa e+ | =, (10) 
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C 2 42,12 
v(t nena) =0, (11) 


where a prime denotes the derivative with respect to r. The matter fields near the boundary 
r — co should behave as 


(2 
g= V+ ais (12) 
aa ae (13) 


where y® = (Ôi, i = 1,2, up to a normalization constant, and Ô; are the corresponding 
dual operators of w in the field theory side. The conformal dimensions of the operators 
are Ay = 24+ V74+m?. u and p are the corresponding chemical potential and charge 
density in the dual field theory, respectively. 

There are only four independent parameters at the tip, i.e., ro, Y(0), 6(70), C (ro). How- 
ever, the above equations of motion have two useful scaling symmetries [26] 


rar, (ee) a, 1) 0% $ > ad, (14) 


C+C-2n6, t>ßt, 4434/8. (15) 


These two symmetries allow us to pick any values for the position of the tip ro and C (ro). 
In our numerical calculations, we will simply choose rp = 1,C (ro) = 0. Furthermore, to 
recover the pure AdS boundary, we also need A(oo) = 0 and C(co) = 0. In general, a 
solution will have a nonvanishing C(co). In this case we will use the scaling symmetry (15) 
to shift C(co) to zero. 


3 Insulator/Superconductor Phase Transition 


From the discussion in Section 2, for given m?,q,w(ro), we can solve the equations of 
motion for the system by choosing (79) as a shooting parameter. After solving the coupled 
equations, we can obtain the condensate (O) from (12) and read off the chemical potential 
u and charge density p through (13). 

In this paper we only focus on the case m? = -5 and q = 2. It is well known that 
in five-dimensional spacetime, when —4 < m? < —3, the scalar field admits two different 
quantizations related by a Legendre transform [27]. Y® can either be identified as a source 
or an expectation value. In this paper we set w“@) = 0 and consider w) acting as the 
vacuum expectation value of the operator Oy. 

It is convenient to make a coordinate transformation from r—coordinate to z—coordinate 
by defining z = 1/r. Therefore, the infinite boundary is now at z = 0 and the soliton tip 
sits at z = 1/ro = 1. Figure(1) exhibits a typical solution, where A(z) and B(z) are 
two metric functions, and (z) and ¢(z) are the scalar field and static electric potential, 
respectively. Note that the metric functions A(z) and B(z) will be used in calculating the 
holographic entanglement entropy. 

For different choices of w(ro), the identification length T in x coordinate will be different. 
In order to compare different solutions with the same boundary behavior, we can use the 
symmetry (14) to set the boundary geometry the same. Taking advantage of the scaling 
symmetry (14), the relevant quantities scale as follows 


A 


1 
Tcr, poop, pap, (0a 0): (16) 
Q 


The combinations pT, pl? and (Oz) 2T are all scale invariant. By using these scale invariant 
quantities, we can accomplish our purpose easily. As in our units setup, the identification 
length IT in the pure soliton is 7, we scale all I for each solution to be I = 7 from now on. 
We should stress here that after the scaling transformation, the tip rọ will be no longer at 
To = 1. 

Figure(2) reproduces the results in Ref. [26] by a slightly different manner. It shows the 
condensate and the charge density as functions of chemical potential. It is clear from the 
figure that as the chemical potential u exceeds a critical value ue for given mass and charge, 
the condensation of the operators emerges. This can be identified as a superconductor 
(superfluid) phase. However, when less than ue, the scalar field is vanishing and this can 
be thought of as the insulator phase, since this system has a mass gap ~ 1/T due to the 
identification in x direction. Note that it is a second order phase transition in our choice 
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Figure 1: A typical soliton solution with nonvanishing scalar hair. Here the value of the 
scalar at tip is wo = 0.5 and the corresponding identification I œ 3.0085. 


of parameters. In addition, let us mention that in the theory described by the action (3), 
except for the AdS soliton solution (insulator) and hairy soliton solution (superconducting 
soliton) mentioned above, there do exist another two solutions: the black hole solution 
without hair (metal phase) and black hole solution with hair (another superconducting 
phase). The complete phase diagram has been drawn in [26]. The details depends on 
the coupling constant q. In this paper we just focus on the case with zero temperature 
insultor/superconductor phase transition. As shown in [26], with our choice of parameters 
m = —2 and q = 2, the phase transition from the insulator to the superconductor always 
happens. In particular, for any values of chemical potential beyond the critical one, the 


superconducting phase always exists. 


4 Holographic Entanglement Entropy 


After finding the gravitational soliton solution with full back reaction of matter fields, 
we are now ready to calculate the holographic entanglement entropy in this holographic 
model. Due to the fact that the choice of the subsystem A is arbitrary, we can define infinite 
entanglement entropies accordingly. In this paper we first consider a simple case where A 
is chosen to be a half of the total boundary space. We assume that the subsystem is defined 
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Figure 2: The condensate of operator (O) (left plot) and the charge density p (right plot) 
versus chemical potential u respectively. Here m? = —15/4, q = 2 and I is scaled to be 
m. The critical chemical potential in this case is ue ~ 0.9442. It is a typical second order 
phase transition which can be seen clearly from u ~ p plot. 


by x > 0 and extended in the y and y directions, where -2 zyx 2 (R> œ), 0<x<T. 
Then the entanglement entropy can be deduced from the formula (2) as 


SA = TER A re > dr. (17) 
For the pure AdS soliton solution, [ = zL and we can read A(r) from metric (4) that 
A(r) = 0. So the entropy can be easily calculated as follows. 
RT fe RT TLR l, LR Rr/1 
grat — dr = ——r? |!" = =Y — —— r = — (= - 1), 18 
A AGa da BG a gG r ra h 


where r = + is the UV cutoff. Note that in our units setup L = 1,r9 = 1, = r. The first 
term is divergent (e — 0) and represents the “area law” [20]. In contrast, the second term 
does not depend on the cutoff and thus is physical important. In fact, this finite term is 
the difference between the entropy in the pure AdS soliton and the one in the pure AdS 
space, because in both cases the divergent part of the entanglement entropy is the same. 
The result in (18) implies that the entropy in the AdS soliton is less than the one in the 
pure AdS space. 

Note that the UV behavior of S4 will not change after the operator condensation. Thus 
the divergent part of the entropy known as the area law will not change because this part 
is only sensitive to UV quantities. This can be understood from the gravity solution side, 
the new solution after the operator condensation still asymptotically approaches to AdS 
space near the AdS boundary. We can write the entanglement entropy in general in the 
half embedding case as 


RT fe aw Rr 
gras _ FO _ fm l 
A = IGu Jn ap 


(19) 


where s has no divergence and s = —1 corresponds to the pure AdS soliton. The value of s 
as a function of chemical potential jz in the superconductor phase is presented in Figure (3). 
One can see from the figure that the entanglement entropy is a constant in the insulator 
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Figure 3: The entanglement entropy as a function of chemical potential in the half em- 
bedding. The solid red curve denotes the entropy in the superconductor phase, while the 
dashed blue line is for the pure AdS soliton solution. The entropy rises as the chemical 
potential u is increased after the phase transition, arrives at its maximum at a certain p, 
and then decreases monotonously. 


phase. After condensate, the entropy first rises and arrives at its maximum as the chemical 
potential u increases, then it decreases monotonously. Thus, from ue to a certain value of 
L, there exits a region where the entropies in the superconductor phase are larger than the 
one in the insulator phase. This behavior of entanglement entropy is quite different from 
the one discussed in Ref. [24], where the case for the holographic metal/superconductor 
phase transition has been studied in detail. It is shown that the entanglement entropy after 
condensate is always lower than the one in the metal phase and decreases monotonously 
as temperature is lowered, and the entropy as a function of temperate is found to have 
a discontinuous slop at the transition temperature T}, which indicates a reduction in the 
number of degrees of freedom after condensate. Similarly, we here also find a discontinuity 
in the slope of the increasing entropy at the critical chemical potential 4e indicated by 
the vertical dotted black line in Figure (3), which may be considered as a significant 
reorganization of the degrees of freedom of the system. Furthermore, this discontinuity 
can be regarded as the signature of the second order phase transition. 

We next consider a more nontrivial geometry which is a straight belt with a finite width 
L along the x direction. The subsystem A sites on the slice r = 4 and expends in y and x 


directions. The holographic dual surface y4 is defined as a three-dimensional surface 
R R 
(=U, 2=c("), =7 <¥ <5 (E>), Us x <P: (20) 
R is the regularized length in y direction. The holographic surface 7,4 starts from z = £ at 


= +, extends into the bulk until it reaches r = r,, then returns back to the AdS boundary 
r=catxr= —£. It is easy to obtain the induced metric on y4 as 


syd 1 dx 
ds? = hida'dx) = (= +r? (=) Jar + r?°dy? + r2e4™ B(r)dx. (21) 


Using the proposal (2), we need to minimize the following functional 


e T ffn “P /1F Br) (dx/dr Pdr. (22) 


We can deduce the equation of motion for the minimal surface from (22) 


Se“ B(r) (dx /d y as 
ree? B(r)(de/dr) _ „3,462 


T+r'B(r)(dcj/drp By: me 


where r, is a constant. 
We are interested in the case that the surface is smooth at r = r,, i.e. dx/dr|,=,, — 00. 
This condition ensures that r, = rą. The width £ of the subsystem A and r, are connected 


by the relation 
z z 1 
=f sar = f — re AG dr. (24) 
2 Ir, dr rori BO e 1) 


r$ B(rx)eAC*) 


Finally, we obtain the entanglement entropy as 


RT f? r4,/B(r)eA _ RT 1 


Sa = — NV S o a ea 
45 n J. JABO O ABN Gr Et) 


(25) 


where the UV cutoff has been taken into consideration. However, there is also a discon- 
nected solution describing two separated surfaces that are located at x = +£, respectively. 
This configuration is just two copies of the half embedding solution that we discussed 
above. So the entanglement entropy of the disconnected configuration is trivially twice of 
the half embedding case (19). 

We find the behaviors of the entanglement entropy as a function of the belt width 
l are quite similar for various chemical potential u or condensation (O). The result for 
u = 2.2016 is shown in Figure (4) as a typical example. When ¢@ > Zmar, there does 
not exist connected surface. The physically favored solution is replaced by the trivial 
disconnected one in this case. On the other hand, when £ < faz, we have three different 
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Figure 4: The entanglement entropy as a function of strip width @ for u = 2.2016. The 
dashed magenta and solid red curves come from the connected solutions, while the solid 
blue one comes from the disconnected one. The lowest curve is physically favored compared 
with others because it has minimal entropy. 


branches of entropy for a given £ , i.e., the upper connected branch (dashed magenta 
curve), the lower connected branch (solid red curve) and the middle disconnected branch 
(solid blue line). The physical entropy is determined by the choice of the lowest one. As 
we can see from Figure (4), for small £, the lower connected branch is favored. However, 
there is a critical value Copiticay above which the lower connected branch is not a physical 
one since its entropy is larger than the one for the disconnected solution. Thus, when £ 
increases, a phase transition occurs as belt width £ = Ceriticgiz. This is just the so-called 
“confinement /deconfinement” phase transition discussed intensively in Refs. [20, 21, 22, 
23]. The transition can be understood as the fact that no correlation could be contributed 
among a distance larger than F due to the mass gap ~ 1/T. 

More numerical results for entanglement entropy are exhibited in Figure (5) with differ- 
ent chemical potential. For 0 < leritical at fixed chemical potential, the entropy is dominated 
by the connected surface and exhibits a non-trivial dependence on £. On the contrary, it 
is governed by the disconnected configuration and is independent of @ for 0 > Coniticat- 
The fact that the phase transition always appears at critical length leritical is observed for 
all values of u matches the prediction in Ref. [21] that this “confinement /deconfinemnet” 
phase transition must be exhibited in any consistent gravity dual of a confining theory. 
We find that max œ (0.220 ~ 0.237)T and leriticat ~ (0.204 ~ 0.221)T for various chemical 
potentials up to u = 5.164 in our numerical calculation. If the belt width is large enough 
(£ > 0.2371), the disconnected branch is always favored for all values of u. Thus, the 
behavior of s as a function of u for fixed £ is the same as the one shown in Figure (3). 


10 


Figure 5: The entanglement entropy as a function of belt width £ for different chemical 
potential u. The dotted black curve is for the pure AdS soliton solution. With a fixed £, 
the entanglement entropy increases (left plot) and then decreases (right plot) as u becomes 
large. 


It is helpful to slice the data in superconducting phase differently. We first study how 
the entropy of the system changes with chemical potential u for fixing belt width £, which 
is presented in the left plot of Figure (6). It can be seen clearly that the entanglement 
entropy as a function of chemical potential for fixed @ has similar behavior as in the half 
embedding case. In particular, the purple curve in the confined phase is the same as that 
in Figure (3), in which £ —> oo. As the chemical potential u increases, the entropy first 
rises and reaches its maximum at a certain value of chemical potential denoted as Umax, 
then it decreases monotonously. The curve will flatten out as we decrease the belt width 
é. The curve will finally become a line in the limit  —> 0, which can be understood from 
the geometric viewpoint. The minimal surface with small belt length can only probe the 
bulk sufficiently near the boundary r — oo where the metric solution is nearly pure AdS. 
Furthermore, the value of Umax also decreases as we decrease the belt width £. The values 
of critical belt width and entropy at the “confinement /deconfinemnet” transition point 
for different chemical potential are drawn in the right plot of Figure (6). The knot is a 
consequence of the non-monotonic behavior of the entanglement entropy. 

In the deconfinement phase with @ < 0.221, the connected branch appears, so the 
quantitative behavior of s(u,/) changes. In order to measure the degrees of freedom at 
energy scale ~ 1/2 in this situation, we would like to calculate the so-called entropic 
c-function first introduced in Ref. [20] 

3 ds(0) 
Coat PTA (26) 
The entropic c-function as a function of belt width @ for a fixed p is plotted in Figure (7). 
Just as it is expected, the value of c decreases monotonically as we increase the belt width £. 
Due to the appearance of mass gap ~ 1/7, the c-function jumps to zero at the critical width 
Lleritical (Ceritical & 0.661 in Figure (7)) and keeps vanishing for larger @. This is just the rea- 
son that one calls the phase transition mentioned above the “confinement /deconfinement” 


11 


—0.8 —0.8 
—1.0F -10 
=1.2+ 
-14}+ 
1 2 3 4 5 å 0.65 0.66 0.67 0.68 0.69 i 


Figure 6: The entanglement entropy of the system changes with chemical potential u for 
fixing belt width £ (left plot). The red curve is for 0 = 0.52 < Cepticay and the blue one 
for £ = 0.57 < Leriticat. The purple curve is for sufficiently large 0 > ln, in the confining 
phase. The right plot presents the behavior of entropy as a function of belt width at critical 
point of the “confinement /deconfinement” transition for different chemical potential. 


phase transition [20, 21, 22, 23]. We also show the change of entropic c-function as we tune 
the value of chemical potential for a fixed small / in Figure (7). We can also see from the 
figure the non-monotonic behavior with the increment of chemical potential. Furthermore, 
the curve becomes flat as we decrease the belt width £. It can be expected that the curve 
will finally become a line in the limit £ —> 0 from the perspective of gravity side. The 
behavior of entropic c-function as a function of chemical potential for a fixed small £ is 
quite qualitatively similar to the entropy in Figure (6). The key difference here is that 
the entropy increases as we increase the belt width, while the c-function decrease for fixed 
chemical potential. But an important feature we observed here is that in both phases, “con- 
finement and deconfinement phases”, the entanglement entropy is always non-monotonic 
in the superconducting phase. 

To summarize, there do totally exist four phases probed by the entanglement entropy 
for a belt geometry for this holographic system, i.e., the confinement /deconfinement in- 
sulator phase, and confinement/deconfinement superconductor phase. These phases are 
characterized by the chemical potential and belt width. In particular, the belt width con- 
trols the “confinement /deconfinement” phase transition. Strictly speaking, the term phase 
transition here is inappropriate since the system itself, i.e., the state of the boundary field 
theory, does not change at all as one changes the belt width @. However, this behavior is 
reminiscent of many thermodynamic phase transitions in holographic calculations (see, for 
example Refs. [28, 29]) and so here we adopt the terminology “phase transition” to convey 
this picture as in the literature [20, 21, 22, 23]. Then one can draw the “phase diagram” as 
in Figure (8) for the system. The left part of the diagram is not new but just the one for 
the pure AdS soliton case which has been studied in Ref. [20]. In the superconductor phase, 
we find a non-trivial phase boundary between the confinement and deconfinement phases. 
In this sense our work is a generalization of Refs. [20, 21, 22, 23] to the superconducting 


12 


C 
a ii 
u=1.1096 ' 
0.4: ! 
0.2 | 
QÜ Sen {| __ 
oe a oe n 2 3 4 5 H 


Figure 7: The entropic c-function as a function of belt width @ for u = 1.1096 (left plot) 
and of chemical potential u at fixed / = 0.5 (right upper curve) and ¢ = 0.6 (right down 
curve), respectively. 


soliton solution. 

The non-monotonic behavior of the entanglement entropy in the superconductor phase 
looks strange at first glance. It is quite different from the case in the metal /superconductor 
transition [24], in which the result of the entanglement entropy decreasing as the tempera- 
ture is lowered agrees with the fact that more degrees of freedom will be condensed in lower 
temperature. In our case, one can see from Figure (6) and Figure (7) that the behavior 
of both entropy and entropic c-function can be affected by belt width. However, the belt 
width £ does not play the essential role for the non-monotonic behavior. As we observed 
above, the non-monotonic behavior occurs in both the confinement and deconfinement 
superconducting phases. 

Due to the lack of the knowledge of the dual field theory in the “bottom-up” approach, 
we do not know the details of the microscopic degrees of freedom and thus it is quite difficult 
to give a definite interpretation of the non-monotonic behavior of the entanglement entropy 
in the holographic insulator/superconductor phase transition. In order to understand the 
numerical result, however, we here try to give a simple physical picture as follows. Note 
that the entanglement entropy as a function of temperature near the critical temperature 
T, shows the same behavior in s-wave and p-wave holographic metal/superconductor mod- 
els [24, 30]. The metal phase can be thought of as the one filled with free charge carriers, 
such as electrons. At the critical temperature the condensate turns on and the free charge 
carriers are continuously condensed to a certain microstate as temperature is lowered. Asa 
result, the entanglement entropy in the superconducting phase is always less than the one 
in the metal phase and decreases monotonically as temperature is lowered. In contrast, 
the insulator phase can be considered as having a little degrees of freedom and almost all 
charge carriers are bound and can not move freely. This is due to the appearance of a mass 
gap ~ 1/T which is the crucial distinction between the holographic metal/superconductor 
and insulator/superconductor transition. At the beginning of the phase transition the 
condensate emerges. Thus, a new kind of charge carriers (some quasi-particles like cooper 
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Figure 8: The “phase diagram” of entanglement entropy for a belt geometry in the holo- 
graphic insulator/superconductor model for m? = -5, q = 2, and [ = z. The phase 
boundary between the confinement phase and deconfinement phase is denoted by the dot- 
ted blue line and solid red curve, while the insulator phase and superconductor phase are 
separated by the vertical dashed line. 


pairs) that can move unconstrained make contribution to the degrees of freedom, which 
leads to the increase of the entanglement entropy. When the chemical potential is enough 
large, however, all bounded carriers are condensed to form cooper pairs. The system in 
this case is full of “cooper pairs” and thus is an ordered state seen from momentum space. 
As a result, the number of degrees of freedom now are even less than the one in the in- 
sulator phase. The above argument may account for the non-monotonic behavior of the 
entanglement entropy in the insulator/superconductor phase transition. 

To further understand the non-monotonic behavior of the entanglement entropy in the 
insulator/superconductor phase transition, it might be helpful to recall the microscopic 
model presented in [9], where the authors wanted to compare the holographic insula- 
tor/superconductor transition with high-T, cuprates, in particular with the RVB (resonat- 
ing valence bond) scenario of high-T, electron doped cuprates (e.g., Ndg_,Ce,CuO,). The 
phase diagrams for both systems are quite similar. The high-T, cuprate at x = 0 (in the 
under-doped region, here x stands for the electron doping) is the antiferromagnetic insula- 
tor. By doping holes, one frustrates the antiferromagnetic order and eventually destroys it. 
Beyond some critical doping £e, the superconducting ground state emerges. The physics 
of the under-doped cuprates is well-described by the t-J model [9]. In the slave-boson 
(SU(2) slave-boson) approach, the electron operators can be decomposed into bosonic and 
fermionc parts, called holons and spinons, respectively. In the confined insulator phase, 
spinons and holons are completely glued together. At the beginning of phase transition, 
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the holons get condensed and the degrees of spinons might be released. This might be 
the reason that the entanglement entropy increases compared to the case in the insula- 
tor phase. When the electron doping continues to increase (chemical potential gets large 
enough in the holographic insulator/superconudctor phase transition), the spinons will 
form cooper pairs and the degrees of freedom associated with spinons get condensed. In 
that case, the entanglement entropy of the system is even less than the case in the insu- 
lator phase. This model might explain the non-monotonic behavior of the entanglement 
entropy. But one caveat should be stressed here that in the gravity setup of the holo- 
graphic insulator/superconductor model, an extra scalar field which is dual to the fermion 
bilinear in the field side does not included [9]. But one expects that including the extra 
scalar field to the model will not change the background geometry very much, so that the 
entanglement entropy keeps the same behavior. Of course, the above discussions are just 
some speculations, further more deep understanding for the non-monotonic behavior of the 
entanglement entropy is called for. 


5 Summary and Discussions 


By using entanglement entropy as a probe, in this paper, we investigated the holographic 
insulator/superconductor phase transition in a five dimensional Einstein- Maxwell-complex 
scalar theory with a negative cosmological constant, we calculated the entanglement en- 
tropy of dual field theory with two kinds of geometry in the AdS boundary: one is a 
half space, and the other is a belt geometry with width £. In the former case, we found 
that there does exist a discontinuity of the slop of the entropy with respect to the chem- 
ical potential at the phase transition point. This discontinuity could be thought of as a 
significant reorganization of the degrees of freedom of the system and a signature of the 
phase transition. Beyond the critical chemical potential, we found that the entanglement 
entropy is not monotonic in the superconductor phase: at the beginning of the transition, 
the entropy increases and arrives at its maximum at a certain chemical potential, and then 
decreases monotonically. 

In the belt geometry case, the so-called “confinement /deconfinement phase transition” 
happens in both the insulator and superconductor phases. As l < ¢¢niticqy the connected sur- 
face is favored and the entropy depends on belt width £ non-trivially, while as 0 > oriticat 
the disconnected surface is favored and the entropy is Z independent. The results en- 
hance the prediction made in Ref. [21] that this phase transition must happen in any 
consistent gravity dual of a confining theory. Thus, as a result, there do exist four differ- 
ent “phases” characterized by chemical potential u and belt width Z, in the holographic 
insulator/superconductor phase transition. The complete “phase diagram” is drawn in 
Figure (8). 

We observed the non-monotonic behavior of the entanglement entropy with respect to 
the chemical potential u in both the “confinement /deconfinement superconductor phases. 
The entanglement entropy firstly rises and reaches its maximum at a certain chemical 
potential and then decreases monotonically (see Figure (3) and Figure (6)). Although 


15 


the entanglement entropy and entropic c-function depend on the belt width, it does not 
play any essential role in the non-monotonic behavior of the entanglement entropy. The 
non-monotonic behavior of the entanglement entropy could be thought of as a key char- 
acteristic of difference between the holographic insulator/superconductor transition and 
the metal/superconductor transition. In the latter case, the entanglement entropy as a 
function of temperature near the critical transition point always monotonically decreases 
as temperature is lowered [24, 30]. Clearly it would be of some interest to check whether 
the non-monotonic behavior of the entanglement entropy is universal in other holographic 
insulator/superconductor model, for example, the p-wave case [32]. 

In addition it is worth stressing here that we did not have a clear microscopic inter- 
pretation for the non-monotonic behavior of the entanglement entropy, but we did present 
some arguments which try to give a physical picture for this behavior. Needless to say, it 
is quite necessary to deeply understand this behavior further. 

In this paper, we have limited ourselves to the choice of parameters as =) = 0, m? = 
—2, and q = 2. Note that the additional complications appear for sufficiently small charge 
q. It is worthy studying the cases with other (m?, q)’s and the case at finite temperature. 
It is also interesting to study entanglement entropy in other holographic superconductor 
models, such as the holographic p-wave superconductor [31, 32, 33, 34]. What’s more, it 
is of some interest to see the effect on the entanglement entropy of the magnetic field in 
the holographic superconductors [35, 36, 37]. We will leave them for future studies. 
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